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Abstract 

The status of relativistic nuclear many-body calculations of nuclear systems to be built up in 
terms of protons and neutrons is reviewed. In detail, relativistic effects on several aspects of 
nuclear matter such as the effective mass, saturation mechanism, and the symmetry energy are 
considered. This review will especially focus on isospin asymmetric issues, since these aspects are 
of high interest in astrophysical and nuclear structure studies. Furthermore, from the experimental 
side these aspects are experiencing an additional boost from a new generation of radioactive beam 
facilities, e.g. the future GSI facility FAIR in Germany or SPIRAL2 at GANIL/France. Finally, 
the prospects of studying finite nuclei in microscopic calculations which are based on realistic 
AW interactions by including relativistic effects in calculations of low momentum interactions are 
discussed. 
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I. INTRODUCTION 



Nuclear systems are very intriguing many-particle systems, which have to be treated by 
means of quantum theory In this article we will consider nuclear systems at low energies to 
be built up in terms of protons and neutrons, ignoring to a large extent all features which are 
related to the substructure of the baryons. As compared to other areas for the application 
of quantum many-body or condensed matter theory, nuclear systems are of finite size. The 
success of the Weizacker mass formula in describing the binding energies of nuclei shows 
that the underlying nucleon-nucleon (NN) forces are of short range and that these binding 
energies result from an interplay between volume and surface contributions. 

The success of this mass formula also allows to complement the nuclear systems of finite 
size, i.e. the nuclei, by a theoretical construct, the system of homogeneous infinite nuclear 
matter. Assuming that there is no Coulomb repulsion between the nucleons one can con- 
struct isospin symmetric nuclear matter assuming an infinite system of constant density with 
the same number of protons and neutrons. The basic properties of this symmetric infinite 
matter can be deduced from finite nuclei by extrapolation. The energy per nucleon should 
correspond to the volume term in the Weizacker mass formula, i.e. the energy per nucleon 
of nuclei with identical number of protons and neutrons, ignoring surface and Coulomb con- 
tributions. This energy per nucleon of about -16 MeV must be obtained for the optimal 
density, i.e. the saturation density po, the density at which the energy per nucleon is mini- 
mal. This saturation density can be deduced from experimental data by extrapolating the 
baryon densities reached in the center of heavy nuclei. This leads to a value of about 0.16 
nucleon/fm -3 . So one first goal of nuclear many-body calculations is to reproduce these 
empirical data of the saturation point for symmetric infinite nuclear matter. 

However, infinite nuclear matter is not only a theoretical construct. The collapse of stars 
after burning all the nuclear fuel in the fusion processes can lead to a nuclear system with 
a homogeneous density over ranges which are so large as compared to the typical distances 
between the nucleons, that surface effects can be ignored. Since, however, these real systems 
of infinite nuclear matter can not switch off the electro-magnetic interaction, the positive 
charges of the protons must be compensated by a corresponding number of electrons or 
muons with negative charge leading to a system in equilibrium with respect to /3-decay, which 
consists predominantly out of neutrons with a small admixture (typically a few percent) of 
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protons. Therefore these objects are called neutron stars. It is one of the challenges of 
nuclear many-body theory to predict the energy density of these infinite nuclear systems 
with large proton-neutron asymmetries over a wide range of densities. The equation of state 
resulting from such calculations is an important ingredient for the simulation of neutron 
stars or astrophysical phenomena like supernovae. 

Of special interest in this context is the study of the crust of such neutron stars. At the 
low densities in the outer crust of neutron stars one can expect the baryonic matter to be 
described as a lattice of "normal" nuclei, whereas the baryonic matter with homogeneous 
density should be formed only at larger densities (larger than about 50 percent of the 
saturation density of symmetric nuclear matter) in the inner crust of the neutron star. 
The transitional region between the isolated nuclei and the homogeneous matter is a very 
interesting example for the phase transition of a quantum- liquid from the droplet phase to the 
homogeneous phase. For this transitional region intriguing structures have been predicted 
like quasi-nuclear bubbles, strings or layers embedded in a sea of neutrons. Because of such 
structures this transitional region has been called the pasta-phase of neutron stars 

These introductory remarks shall demonstrate that the nuclear many-body problem is 
indeed a very interesting topic in theoretical physics. It is of interest on its own, because we 
want to understand the bulk properties of nuclei in terms of the underlying NN interaction. 
Of special interest is the study of exotic nuclei with a large neutron excess and the extrapo- 
lation to the matter in neutron star, since this provides the information for the simulation of 
astrophysical processes. Theoretical studies of the nuclear many-body problem in terms of 
nucleonic degrees of freedom are also of large interest because they may finally exhibit the 
limitations of such nuclear structure studies ignoring the sub-nucleonic degrees of freedom 
even in describing low energy properties of nuclei. 

Therefore a lot of effort has been made during the past decades to solve this conventional 
nuclear many-body problem. In general one may distinguish between so-called microscopic 
or ab initio approaches on one hand and more phenomenological approaches on the other 
hand. 

In what we will call the phenomenological approach, one starts with a simple model for a 
NN or even 3N interaction and adjusts the parameter of this model to describe the saturation 
point of symmetric nuclear matter and properties of a few selected nuclei. A typical example 
of such a phenomenological interaction is the family of Skyrme interactions 4H6|. In this 
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way one can describe rather accurately properties of nuclei in detail. However, there may be 
a problem in the predictions derived from these phenomenological studies for the systems 
at large densities and proton-neutron asymmetries, which are relevant for the study of the 
astrophysical objects. 

In the microscopic approach the parameters of the NN interaction are adjusted to de- 
scribe the experimental data of NN scattering including energies up to the threshold for 
pion production and the properties of the deuteron. The NN scattering data, in particular 
the change of sign in the S-wave phase shifts suggest that such a realistic NN interaction 
should contain repulsive components of short range and attractive components with medium 
and long ranges. In fact, some of the early models for a local NN interaction considered a 
hard core potential, with infinite repulsion for relative distances below r c w 0.4 fm {?], Isj]. 
Also modern versions for such a local NN interaction contain a considerable amount of 
cancellation between very repulsive components of short range and attractive ones at larger 
ranges {9]. 

But also other models, like the one-boson-exchange model 10], which is based on the 
meson exchange theory for the NN interaction and leads to a non-local potential, contain 
rather strong components which are compensating each other to a large extent. 

Employing such realistic NN interaction with very strong components it turns out, that 
perturbative approaches in terms of the NN interaction do not lead to reliable results. In 
fact, using such interactions within a mean field or Hartree-Fock approach does not provide 
binding energy for nuclear matter Therefore attempts have been made to sum up 

all ladder diagrams and evaluate the energy of nuclear matter in terms of the T matrix. 
Hartree-Fock calculations in terms of this scattering matrix lead to a sufficient amount of 
binding energy for symmetric nuclear matter but fail to produce a saturation point: the 
energy per nucleon calculated in this way decreases with increasing density and does not 
provide a minimum 10]. 

The situation is improved if one calculates the ladder diagrams accounting for nuclear 
medium effects in the sense that the intermediate two-nucleon states are restricted to the 
states above the Fermi momentum and by accounting for the binding effects in two-nucleon 
propagator. This leads to the so-called Brueckner G-matrix and the Brueckner-Hartree-Fock 
(BHF) approximation. 

The Pauli- and dispersive features included in the G-matrix lead to a minimum in the 
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calculated energy versus density plot, which means that they provide a theoretical result for 
the saturation point of nuclear matter. A lot of effort has been made to generate a realistic 
NN interaction, which using the BHF approximation also reproduces the experimental result 
for the saturation of nuclear matter. However, all such calculations produced saturation 

n 

points on the so-called Coester band [12]. NN interaction which are rather soft tend to 
produce energies per nucleon which are close or even below the empirical value of -16 MeV 
but at a saturation density which about twice or even larger than the empirical po- Other 
interactions, which are stiffer and contain larger tensor components may predict the correct 
saturation density but yield a value for the energy per nucleon, which is considerably large 
than the empirical value (-10 MeV compared to -16 MeV). Anyway none of the interactions 
predicts a saturation point, which is close to the experimental result. 

Attempts have been made to go beyond this BHF approximation. It is one drawback of 
this approach that particle-particle and hole-hole ladder terms are treated differently. This 
implies that the BHF approximation is not symmetry conserving in the sense of the Kadanoff 
and Baym 13[], which implies e.g. that it does not fulfill the Hugenholtz van Hove theorem. 
This has been improved by completing the particle-particle ladders by corresponding hole- 
hole terms within the ladder approximation of the self-consistent Greens function (SCGF) 
approach [ijj]. It turned out however, that the inclusion of the hole-hole terms does not 
lead to a significant change for the predicted saturation properties of nuclear matter. This 
supports the general strategy of the Brueckner Bethe hole line expansion. 

The next step in the Brueckner hole-line expansion beyond the BHF approximation, is to 
account for all three-particle ladders by solving the Bethe- Fadeev equation 15] . Techniques 
and details, how to solve the Bethe- Fadeev equation have been described by Day nu- 
merical calculations for including the three-hole line terms have been presented by Song et 
al. [l7]. These studies show that the inclusion of three hole line terms provides the nice 
feature that the results are getting essentially independent on the choice for the particle 
spectrum, a choice which can not be decided within the the two hole line approach. How- 
ever, also the three hole line terms yield only a small correction as compared to the BHF 
approximation and therefore the predicted saturation points still lead to a Coester band, 
which fails to come close to the empirical result. 

A similar feature can also be observed in the calculation of bulk properties of finite nuclei 
using the BHF approximation applied to realistic NN interactions [ill ]. Such calculations 
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tend to predict to little binding energy and radii for the charge distribution, which are to 
small as compared to the experimental findings. 

The same is true for other many-body techniques, which have been used to evaluate 
the properties of infinite nuclear matter or finite nuclei employing realistic iViV interac- 
tions within a non-relativistic framework. Such many-body techniques include the coupled 



functions 



cluster or exponential S method [18|, |l9j , variational methods allowing for correlated basis 



201 ] and Quantum Monte Carlo Calculations |2l|. Results close to the experimental 



data could only be obtained after introducing empirical three nucleon forces. 

An alternative method, which is also based on realistic NN interactions, has received 
a lot of attention during the last few years: the use of low-momentum interaction V\ ow _^. 
The basic idea of V[ ow _k is to separate the long range or low-momentum components of 
the NN interaction, which are constrained by the NN scattering matrix below the pion 
threshold, from the high-momentum components, which may strongly depend on the un- 
derlying model of the NN interaction. By introducing a cutoff A in momentum space, one 



separates the Hilbert space into a 



malization technique (see, e.g. 
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ow- momentum and a high- momentum part. The renor- 



26]) determines an effective Hamiltonian, which treated 



within the low-momentum space yields the same results than the original Hamiltonian for 
the corresponding states. This renormalization can be achieved by means of the unitary 



model operator approach 
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241 ] . It is a nice feature of this V[ ow _k interaction approach, 



that using a cutoff A of 2 fm -1 , which corresponds to the pion threshold in NN scattering, 
one obtains a renormalized interaction, which is essentially independent on the underlying 
realistic interaction, which was used to generate Viow-k- Therefore one has got a unique NN 
interaction model. 

Another advantage of using this V[ ow _k approach is that the strong short range components, 
which are characteristic of a realistic NN interaction have been integrated out and a Hartree- 
Fock (HF) calculation is a reliable approximation for the solution of the many-body problem. 
The price one has to pay however, is that such HF calculations using Vjow-k interaction 



do not 
matter 



ead to a saturation point when calculating the binding energy of infinite nuclear 
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281 ] . This is very similar to the results obtained in HF calculations using the NN 



T matrix already discussed above [10j. Indeed Vj ow _k can be considered to be rather close 
to the T matrix with the only difference that the NN states below the cutoff are excluded 
from the ladder terms in V[ ow _k as compared to T. Also in this case a possible way out of 
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the saturation problem is to include a phenomenological three-nucleon force. 

All this demonstrates that three-nucleon forces seem to be necessary to obtain the empir- 
ical saturation property of nuclear matter within a non-relativistic theory based on realistic 
models of the NN interaction. 

On the other hand, however, relativistic calculations have been successful to derive the 
empirical saturation point of symmetric nuclear matter from realistic meson exchange mode 



of the NN interaction without the necessity to adjust appropriate three-nucleon forces. 
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411 ] At first sight relativistic effects seem to be negligible when dealing with the nuclear 



many-body problem, as the binding energies of nucleons are very small as compared to the 
nucleon rest mass. The meson exchange model of the NN interaction, however, predicts 
a strong cancellation of various components in the relativistic self-energy of nucleons in 
a nuclear medium. A strong attractive part (around -300 MeV), which transforms like a 
scalar under a Lorentz transformation and mainly originates from the exchange of the scalar 
a meson in realistic One-Boson-Exchange (OBE) models of the NN interaction, is partly 
compensated by a repulsive component, which transforms like a time-like component of a 
Lorentz vector and reflects the repulsion due to the u meson exchange. The effects of this 
two contributions cancel each other to a large extent in calculating the single-particle energy 
of the nucleon, which is indeed small (-40 MeV) on the scale of the mass of the nucleon. 

Inserting the nucleon self-energy in a Dirac equation for the nucleon, the strong scalar 
part of the self-energy yields an enhancement of the small component in the Dirac spinor 
of the nucleon in the medium as compared to the vacuum. This modification of the Dirac 
spinors in the nuclear medium leads to modified matrix elements for the OBE interaction. It 
is this density dependence of the nucleon Dirac spinors and the resulting medium dependence 
of the NN interaction, which moves the saturation point calculated for nuclear matter off 
the Coester band in such a way, that the empirical values for the energy per nucleon and the 
saturation density can be reproduced without the necessity to adjust any parameter or to 
introduce many-nucleon forces. Reasonable saturation properties have even been obtained 
within the framework of the V[ ow _k approach, if this change of the nucleon Dirac spinors in the 
nuclear medium is taken into account [42J. Therefor one may argue that the three-nucleon 
forces, which are required to achieve saturation within a non-relativistic treatment, is an 
attempt to simulate the effects of these relativistic features. 

Although the relativistic many-body approach, the Dirac-Brueckner-Hartree-Fock ap- 
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proach in particular, has been very successful in describing the saturation property of nuclear 
matter, only a few attempts have been made to extend the DBHF approximation also to the 
evaluation of finite nuclei 
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46| . The consistent treatment of correlation and relativistic 



effects for finite systems is a rather involved problem. 

On the other hand, very extensive studies of finite nuclei have been made within relativis- 



tic mean field approximation 
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48j using a parametrization which is adjusted to describe 



the experimental data of finite nuclei. With the advance of the renormalization techniques 
leading to effective low-momentum interactions V[ ow _k the time may be appropriate to con- 
sider the possible bridge between ab initio relativistic calculations based on NN interactions 
fitted to the scattering data and these phenomenological relativistic mean field calculations. 

We hope that this review on the status of relativistic nuclear many-body calculations will 
help to built this bridge. 

II. RELATIVISTIC MANY-BODY APPROACHES 



During the past decades, much success has been achieved in nuclear physics by relativistic 
many-body approaches. These approaches are the Dirac-Brueckner-Hartree-Fock (DBHF) 
approaches and relativistic density functional theories. The DBHF approach is an ab initio 
approach, whereas relativistic density functional theory is a phenomenological approach. 



These relativistic Brueckner ca 
various groups 
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32 



37 
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culations are not straightforward and the approaches of 
4l| are similar but differ in detail, depending on solution 
techniques and the particular approximations made. Therefore, a general description of the 
relativistic Brueckner approach will be given and the implications of the several approxima- 
tion schemes will be discussed. 

One of the most outstanding relativistic density functional theories is the relativistic 
Hartree approach with the no-sea approximation, namely the relativistic mean-field (RMF) 
theory. Therefore, the DBHF approach is treated in Sec. Ill Al the RMF theory will be 
discussed in Sec. Ill B\ and the extensions to the RMF theory in Sec. Ill CI 
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A. DBHF approach 



The main framework of the relativistic DBHF approach consists of a set of coupled 
equations, in which the interaction of the nucleons, the propagation of the particles, and 
self-energies occur. Firstly, the in-medium interaction of the nucleons is treated in the ladder 
approximation of the relativistic Bethe-Salpeter (BS) equation 

T = V + i J VQGGT, (1) 

where T denotes the T matrix, V the bare nucleon-nucleon, and Q the Pauli operator, 
which prevents scattering to occupied states. Furthermore, the Green's function G, which 
describes the propagation of the dressed intermediate off-shell nucleons in Eq. JT]), fulfills 
the Dyson equation 

G = Go + G EG, (2) 

where Go denotes the free nucleon propagator and E represents the self-energy. This self- 
energy E in the Hartree-Fock approximation is given by 

E = -ij (Tr[GT] - GT). (3) 

F 

Therefore, this coupled set of equations, Eqs. ([I])-®, represents a self-consistency problem 
and has to be iterated until convergence is reached. 

The structure of the self-energy follows from translational and rotational invariance, her- 
miticity, parity conservation, and time reversal invariance. Therefore, the most general form 
of the Lorentz structure of the self-energy is given by 

E = E s - 7 A1 E M (4) 

with 

E^ = E u" + E V A^, (5) 

where the E s , E D , and E v components are Lorentz scalar functions, which depend on the 
Lorentz invariants k 2 , k ■ j and j 2 , with j M the baryon current and k^ the nucleon four- 
momentum. Therefore, these Lorentz invariants can be expressed in terms of k , |k| and the 
Fermi momentum kp. Furthermore, the projector in Eq. (jSJ) is given by = gV—u^u" 



with streaming velocity = j M / In the nuclear matter rest frame, = the self- 
energy then has the simple form 

k F ) = T, s (k, fcp) - 70 S Q (fc, fc F ) +7 • k£ v (£;, fc F ). (6) 

These E s , £ Q , and S v components of the self-energy can be easily calculated by taking the 



respective traces 
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42| 



Ss = ^r[S], E o = -Ur[ T0 £], S v = — -Lr[ 7 -k£]. (7) 



The masses and momenta of the nucleons inside nuclear matter are modified by the 
nuclear medium 

m*{k,k F ) = M + ^ei: s {k,k F ), fc*'* = F + 9ffe£"(A;,Jfe F ). (8) 

Written in terms of these effective quantities, the Dirac equation has the form 

\Ji*--m* -i 3m £]«(£;) = 0. (9) 

The imaginary contribution to the self-energy is due to the possible decay of particle states 
above Fermi sea into hole status within the Fermi sea. To simplify the self-consistency scheme 
this decay possibility is neglected and one works in the quasi-particle approximation, i.e. 
3mE = 0. Furthermore, introducing the reduced effective mass, 

m*(k, k F ) = m*{k, k F ) / (1 + £ v (Jfe, k F )) , (10) 

and the reduced kinetic momentum, 

A;; = ^/(1 + £ v (Mf)), (11) 

the Dirac equation written in terms of these reduced effective masses and momenta has the 
form 

[ 7m F m - m*{k, k F )}u{k, k F ) = 0. (12) 
The solution of the Dirac equation in Eq. (fl2l) is 



<<^-M = a/^^ ,1 I Va- (13) 



2m* 2A l k l 

X £*(k)+mj 
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where E*(k) = y k + fh*p denotes the reduced effective energy and xa a two- component 
Pauli spinor with A = ±|. The normalization of the Dirac spinor is thereby chosen as 
u\(k, kp)u\(k, /cf) = 1- From the Dirac equation in Eq. (fl2l) one derives the single-particle 
potential U = 7°S. It can be obtained by calculating the expectation value of U. Therefore, 
one sandwiches U between the effective spinor basis ([TBI . 

v ; (u(k)\u{k)) E*(k) V V ;I 1 V ;/ E*(k) K ' 

However, the results depend strongly on approximation schemes and techniques used in 
the DBHF approach. In principal one can distinguish between two frequently used methods, 
the fit method versus the projection technique method. The fit method is the more simple 
method; one avoids cumbersome projection techniques which are required using the trace 



formulas in Eq. (j7j). This method was originally proposed by Brockmann and Machleidt [29 1 . 
In this approach, the scalar and vector self-energy components are directly extracted from 
the single particle potential U. A fit to the single-particle potential U in Eq. (1T41) delivers 
the density dependent components S s and S . An attempt, which tries to extend this 



method and to extract momentum dependent fields by fitting procedures [50j], suffers from 
large uncertainties, since one then tries to extract two functions out of one. Therefore, only 
mean values for the self-energy components, where the explicit momentum dependence has 
already been averaged out, are relatively reliably obtained. However, the extrapolation to 
asymmetric matter already makes even this procedure ambiguous for isospin asymmetric 
nuclear matter [38], since it introduces two new parameters in order to fix the isovector 
dependencies of the self-energy components. 

The other method, the projection technique method, is rather complicated, but accurate. 



For example, they have been used in Refs. 32 
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40 



41] . It requires the knowledge 



of the Lorentz structure of the positive-energy-projected in-medium on-shell T matrix. The 
T matrix has to be projected onto covariant amplitudes. Therefore, the scalar and vector 
components of the self-energies can directly be determined from the projection onto Lorentz 
invariant amplitudes. Ambiguities [51] arise due to the restriction to positive energy states, 
since pseudoscalar (ps) and pseudovector (pv) components can not uniquely be disentan- 
gled for on-shell scattering. However, these ambiguities can be minimized by separating the 
leading order, i.e. the single-meson exchange, from the full T matrix. Therefore, the con- 
tributions due to the single- 7r and-// exchange are given in the complete pv representation. 
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For the remaining part of the T matrix, the ps representation is chosen 40|, [41 ] . 



B. Relativistic mean-field theory 

A Lagrangian density of interacting many-particle system consisting of nucleons and 
mesons is the starting point of a relativistic mean-field (RMF) theory. The different RMF 
models differ from each other in the mesons and in the various couplings to the nucleon 
field included in the Lagrangian density. The most simple RMF model is the linear a — uj 
model, which only includes linear coupling terms of the cr-meson and of the w-meson. These 
models provide a relatively poor description of nuclear systems. Therefore, nonlinear terms 
are added to the Lagrangian density to improve the description of isospin symmetric nuclear 
systems [52J and isovector mesons are included in the models for the description of isospin 



asymmetric nuclear systems 53|, [54 1 . 



The Lagrangian density of the RMF theory presented here is taken as general as possible 
and includes as well the isoscalar mesons a and u as the isovector mesons 5 and p. The 
electromagnetic interaction is mediated by the photon 7. Therefore, the Lagrangian density 
consists of four parts: the free baryon Lagrangian density the free meson Lagrangian 
density Cm, the interaction Lagrangian density C int and the nonlinear meson Lagrangian 
density C n \. 

£ = £b + Cm + Ant + Cni, (15) 
which takes on the explicit form 



£m=\Y. (d^d^-m 2 ^) 

L=cr,S 



2 



Ant = - 9a^^a^ ~ g S ^T^ ~ &,*7^(L)* 



(16) 



-eki(l + r 3 )AfJ, 



M2V X ^ (7) 
I 1, 

nl 



with the field strength tensor F^ K ^ V = d^Ar^ — d v A^ K ^ for the vector mesons. In the 
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above Lagrangian density the nucleon field is denoted by \1/ and the nucleon rest mass by 
M. The scalar meson fields are $ CT and and the vector meson fields are Ai w \ and A( p y 
Furthermore, the bold symbols denote vectors in the isospin space acting between the two 
species of nucleons. The mesons have rest masses m a , m^, ms, and m p , and couple to the 
nucleons with the strength of the coupling constants g a , g$, g^ , g p , and f p . The coupling 
constants of the first and second order nonlinear cr-meson self-interactions are a and b, 
respectively. Furthermore, the electromagnetic field couples to the nucleons through 
the electron charge e 2 = Ana where a is the fine structure constant. 

To obtain the field equations, we minimize the action for variations of the fields <ft K 
included in the Lagrangian density (eq. [TS]) 

8 jT dtj d z xC{<p K {x),d p <p K {x),t) =0, (17) 

where <p K stands for the nucleon field, the meson fields and the electromagnetic field. Finally 
the Euler-Lagrange field equations are obtained for each field 4> K 

d dC _^^L_ Q (18) 

Evaluating the Euler-Lagrange field equations for each field <p K , we obtain a Dirac equation 
for the nucleons and Klein-Gordon and Proca equations for the meson fields. This Dirac 
equation for the nucleon field can be written as 

(17^ - M - S) * = (19) 

with the self-energy £ generated by the interaction terms 



L (p) 



(20) 



Secondly, the equations for meson fields are given by 

(□ + m 2 )$ a + a$ CT 2 + 6$ CT 3 = -g a W, (21) 

{n + m 2 s )<f> 5 = -98^, (22) 

(□ + ml )A (U)I1 = gjlifl, (23) 

(□ + m 2 p )A {p)p = g^r^ + d"^ttT<v*, (24) 

UA {1)P = e*|(l + r 3 )7 A4 *- (25) 

Although these RMF models provide a reasonable description of nuclear matter and finite 
nuclei, medium effects are still not satisfactorily treated. 
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C. Extensions of the RMF theory. 



In this section, some extensions to the standard RMF theory will be discussed to obtain 
a more satisfactory treatment of the medium effects. Therefore, RMF models were intro- 
duced, which allow for a density dependence of meson-baryon vertices g K and f K with k a 
meson. Due to this density dependence of the coupling functions, the nonlinear terms are 
not needed anymore and are omitted in these theories. The density p(^, \t r ) is obtained from 
the nucleon field This introduced density dependence can improve the capability of the 
model significantly, in particular when it is based on microscopic many-body calculations. 
In literature, a scalar density dependence or a vector density dependence can be found. It 
turns out that the dependence on the zero component of the vector density, the baryon 
density p = Tr(\I/7o\l/), is the most suitable one since it describes finite nuclei better and has 
a natural connection to the vertices in microscopic many-body approaches like the DBHF 



calculations 
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56] 



These density dependent coupling functions can be determined by a fit either to experi- 
mental data or to the self-energies of the Dirac-Brueckner-Hartree-Fock approach for asym- 
metric nuclear matter. To be consistent with more involved theories for nuclear physics 



the second method is pre : 



tic mean-field model 



48 



erred, since in contrast to the purely phenomenological relativis- 



571 ] . the density dependence of the interaction has a microscopic 
foundation when the density dependent RMF (DDRMF) approach is based on microscopic 
many-body calculations. This should give one confidence when the model is used for extreme 
cases such as neutron-rich nuclei. Therefore, the relation between the the DDRMF theory 
and DBHF calculations will be discussed in Sec. E] 

As a consequence of this density dependence of meson-baryon vertices, we have to vary 
the Lagrangian density by 

= d£ dCSp_ 
5^ dp <ftf' 1 ' 

where the second expression generates the so-called rearrangement contribution E^ to the 
self-energies of the nucleon field. This rearrangement contribution YP'^> has to be added to 
the Dirac equation (fTU|) 

(i-y^d" -M-E)^ = 0^ (i 7(U <9" - M - (E + E (r) 7o )) * = 0, (27) 



14 



where the rearrangement self-energy contribution reads 

These rearrangement contributions should be taken into account and are essential to provide 
a symmetry conserving approach, which implies that energy-momentum conservation and 
thermodynamic consistency like the Hugenholtz - van Hove theorem are satisfied 
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Note that such rearrangement terms are not present in the microscopic DBHF approx- 
imation. That is why the DBHF as well as its not-relativistic counterpart, the Brueckner- 
Hartree-Fock (BHF) approximation, are not symmetry conserving. An extension of BHF, 
which accounts for particle-particle as well as hole-hole ladders, would be required to achieve 
this consistency [14| . This implies that the definition of the nucleon-selfenergy in ([3]) must 
be extended by terms of higher order in the hole-line expansion to obtain such a number 



conserving approach 59] 



Another extension of the RMF theory is the inclusion of Fock terms to obtain a relativistic 
Hartree-Fock (RHF) theory. In this approach, an additional meson, the 7r-meson, is usually 
included. This means that the Langrangian density has to be extended, 

Ar = l(d^d^ n - ml$l) - Jz-Vtw^Q*]*, (29) 

and consequently all equations derived from it. The presence of a it field is a typical feature 
of the RHF theory compared to the RMF theory, since this 7r field only contributes to the 
Fock terms. Therefore, it is omitted in the RMF theory. 

III. EFFECTIVE MASS 

The introduction of an effective mass is a common concept in nuclear physics. It is used to 
characterize the quasi-particle properties of a particle inside a strongly interacting medium. 
Furthermore, it is a well-known fact that the effective nucleon mass in nuclear matter or 
finite nuclei deviates substantially from its vacuum value 60M62| . However, the phrase of 



an effective nucleon mass has been used to denote different quantities: the nonrelativistic 
effective mass m* NR and the relativistic Dirac mass m* D . Although these quantities are 
related, they are based on completely different physical concepts. 
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FIG. 1: The Dirac mass (dotted line), nonrelativistic masss rn* NR (dashed-dotted lines), /c-mass 
mt(k) (solid lines), .E-mass m* E (k) (dashed lines) in isospin symmetric nuclear matter as a function 
of the momentum k = |k| obtained from DBHF calculations at a fixed nuclear density of p = 
0.181 fm~ 3 . 



The Dirac mass is a genuine relativistic quantity and can only be determined from rel- 
ativistic approaches. It is defined through the scalar part of the nucleon self-energy in the 
Dirac field equation which is absorbed into the effective mass 

m* D (k,k F ) = M + m2 s {k,k F ), (30) 

where S s is the scalar part of the nucleon self-energy. The Dirac mass accounts for medium 
effects through the scalar part of the self-en ergy and is a smooth function of the momentum 



as shown in Fig. [I] for a DBHF calculation 63]. 



On the other hand, the nonrelativistic mass m* NR parameterizes the momentum depen- 
dence of the single particle potential. Therefore, it can be determined from both, as well 
relativistic as nonrelativistic approaches. The nonrelativistic mass is defined as 

m* NR =\k\[dE/d\k\}~\ (31) 
16 



where E is the quasi-particle's energy and k its momentum. When evaluated at k — kp, 
Eq. fl3TT) yields the Landau mass. In the quasi-particle approximation, i.e. the zero width 
limit of the in-medium spectral function, the energy E of the quasi-particle and its momen- 
tum k are connected by the dispersion relation 



E 



2M 



(32) 



Hence, the equations ( 13T1) and ( )32l) yield the following expression for the nonrelativistic 
effective mass 

m * NR= J4 + W\W\ m{lklU = E) ' ^ 
Since the nonrelativistic effective mass parameterizes the momentum dependence of the 

single particle potential, it also is a measure of the nonlocality of the single particle potential 

U. 

These nonlocalities of U can be due to nonlocalities in space or in time. The spatial 
nonlocalities result in a momentum dependence and are characterized by the so-called k- 
mass, which is obtained from the derivative of U with respect to the momentum at fixed 
energy, 



m k {k) 
M 



1 + 



MdU(k,u) 



(34) 



64 



651 ] and 



k dk 

These spatial nonlocalities of U are mainly generated by exchange Fock terms 
the resultin g k- mass is a smooth function of the momentum as shown in Fig. [T]for a DBHF 
calculation [63]. 

On the other hand, nonlocalities in time result in an energy dependence and are expressed 
in terms of the E-mass, which is given by the derivative of U with respect to the energy at 
fixed momentum, 



m E {u) 
M 



dU(k } u) 
dui 



(35) 



The nonlocalities in time are generated e.g. by Brueckner ladder correlations due to the 
scattering to intermediate off-shell states. These correlations mainly have a short-range 
character and generate a strong momentum dependence with a characteristic enhancement 



of the E-mass slightly above the Fermi surface 



62 



6J-|66]. 



The effective nonrelativistic mass defined by Eqs. (|3"T|) and ( 133]) is given by the product 



of the k-mass and the E-mass 



66] 



m*(k) m *k{k) m* E (uj = E) 



M 



M 



M 



(36) 
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FIG. 2: The effective nonrelativistic (Landau) and Dirac mass in isospin symmetric nuclear matter 
at k = |k| = kp as a function of the Fermi momentum kp for various relativistic models are plotted. 
In addition, the Landau mass of a nonrelativistic BHF calculation is shown. 



Therefore, it contains nonlocalities as well in space as in time and should show a typical peak 
structure around kp. This narrow enhancement near the Fermi surface reflects - as a model 
independent result - the increase of the level density due to the vanishing imagi nary part of 



the optical potential at fcp, which for example is seen in shell model calculations |61l. I62L 166] . 
In order to reproduce this behavior, one should, however, account for correlations beyond 
mean-field or Hartree-Fock. 

The enhancement of the effective mass m* due to the effective E'-mass in Eq. (l36j) is not 
strong enough to compensate for the effects of the fc-mass. Therefore, the final effective 
mass is below the bare mass M as can be seen in Fig. [TJ 



18 



Another aspect is the density dependence of these two different masses. The Dirac mass 
of relativistic models decreases continously with increasing Fermi momentum Jcj _as shown 
in Fig. HJ The plotted results in Fig. [2] are from a relativistic DBHF approach 



67 



68 



based 



on the Bonn A interaction and from RMF calculations of the NL3 and DD model. The 
NL3 model is a nonlinear parameterization 69J that is widely used in RMF calculations. It 
contains the cr-meson, the w-meson, the p-meson without its tensor coupling, the electro- 
magnetic field, and nonlinear self-interactions of the a field. Compared to the Lagrangian 
density of the NL3 model, the DD model is based on a Lagrangian density without nonlinear 
self-interactions of the a field, however, with density dependent meson-nucleon coupling ver- 
tices (70). Initially, the Landau mass also decreases with increasing Fermi momentum k-p like 
the Dirac mass. However, it starts to rise again at high values of the Fermi momentum hp. 
In addition to the results from relativistic approaches, the Landau mass of a nonrelativistic 
BHF calculation, which is also based on the Bonn A interaction, is plotted in Fig. |2j The 
agreement between the calculated Landau mass from the nonrelativistic and the relativistic 
Brueckner approach is quite good. Therefore, the often discussed difference between effec- 
tive masses obtained in the various approaches is mainly due to different definitions, i.e. 
nonrelativistic mass versus Dirac mass. Although the RMF models qualitatively show the 
same behavior as the Brueckner approaches, their masses are lower compared to the ones in 
the Brueckner approaches. 

Another interesting issue is the proton-neutron mass splitting in isospin asymmetric nu- 
clear matter. This issue will be of relevance in the study of drip-line nuclei. Furthermore, 
it can be expected to have important effects on transport properties, like collective flows, of 
dense isospin asymmetric nuclear matter that will be reached in heavy ion experiments. 

The DBHF calculations based on projection techniques predict a Dirac mass splitting of 



m* D n < m* D in neutron-rich nuclear matter 



36 



38 



40 



41 



67 



as can be observed in 



Fig. |3]for such a DBHF calculation [41] . In contrast, the DBHF calculations based on the fit 
method yield a mass splitting of m* D n > m* D p 38, 39, [zi] as shown in the left panel of Fig. [3] 
for such a calculation 39j. In the fit method, momentum independent self-energy compo- 
nents are obtained, since the explicit momentum dependence has already been averaged 
out. Hence, only mean values for the self-energy components are obtained. This method 
is relatively reliable in isospin symmetric nuclear matter. However, the extrapolation to 
isospin asymmetric matter introduces two new parameters in order to fix the isovector de- 
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FIG. 3: Left panel: Results obtained from different DBHF calculations for the Dirac neutron and 
proton effective masses as a function of the density p for a value of the asymmetry parameter 
j3 = pn p Pp = 1. One of the DBHF calculations is based on projection techniques and the other 
calculation is based on the fit method. Right panel: The Dirac neutron and proton effective masses 
of a RMF theory, which includes the isovector 5-meson, are shown. 



pendencies of the self-energy components. This introduction of two new parameters makes 



the fit procedure ambiguous in isospin asymmetric nuclear matter 



38. 



72]. Furthermore, 



RMF theories with the isovector p- and 5-mesons included predict a Dirac mass splitting of 



rrr Dn < m* D [53|, |73| as you can see in the right panel of Fig. [3] for such a RMF calcula- 
tion [53| . This behavior of the neutron-proton mass splitting is in agreement with results 
from the DBHF calculations based on projection techniques. When only the p-meson is 
included, the RMF theory predicts equal masses, m* Dn = m* D . Hence, the 5-meson is 
responsible for the mass splitting in the RMF theory. 

For completeness, attention will now be given to the nonrelativistic effective mass, which 
is obtained in Fig. H] from a BHF calculation (t4 ] by taking the product of the effective 
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FIG. 4: Effective fc-mass m£(fc) (solid lines), effective .E-mass m* E (k) (dashed lines) and effective 
nonrelativistic masss rn* NR (dashed-dotted lines) for neutrons and protons (lines with symbol) as 
obtained from the BHF calculations for a value of the asymmetry parameter f3 = 0.5 at fixed 
nuclear density p = 0.17 fm~ 3 . The Fermi momenta for protons and neutrons are indicated by 
vertical dotted lines. 

/c-mass and the effective E-mass as done in Eq. fl36l) . The effective fc-mass is larger for the 
neutrons than for the protons, whereas the values for the effective E'-mass are larger for the 
protons than for the neutrons in neutron-rich nuclear matter. Since the difference between 
the effective E-mass of the neutron and the proton is not large enough to compensate for the 
isospin asymmetric nuclear matter effects of the fc-mass, the nonrelativistic proton mass is 
in general smaller than the one of the neutron, in particular at k = |k| = kp. This result for 
the proton-neutron mass splitting is a general feature of BHF calculations 
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FIG. 5: Neutron and proton effective mass obtained from a DBHF calculation based on projection 
techniques are plotted as a function of the momentum k = |k| for a value of the asymmetry 
parameter (3 = 1 at fixed nuclear density p = 0.166 fm~ 3 . In addition, the effective mass in 
symmetric nuclear matter is given and the neutron Fermi momenta are indicated by vertical dotted 
lines for isospin symmetric nuclear matter and for pure neutron matter. 



In contrast, the RMF theorie and the DBHF calculations based on on projection techniques 
predict an opposite behavior for the Dirac mass spliting. 

As already mentioned before, the nonrelativistic mass m* NR can be determined from both, 
as well relativistic as nonrelativistic approaches. The nonrelativistic mass derived from the 
RMF theory shows the same behavior as its Dirac mass, namely m* NRn < rn* NRp [l^. The 
nonrelativistic mass derived from the DBHF approach based on projection techniques shows 
a nonrelativistic mass splitting of m* NRn > m* NRp except around the peak slightly above 
the proton Fermi momentum k-pp- This behavior is opposite to its Dirac mass splitting 
of m* D n < m* n , as can be seen in Fig. [5] for a DBHF calculation based on projection 



techniques 40|, |67|, |68) . However, it is in agreement with results from nonrelativistic BHF 
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calculations [64|, |65|, |75|. Since these masses are based on completely different physical 
concepts, the opposite behavior between the Dirac mass splitting and the nonrelativistic 
mass splitting is not surprising. On the one hand, the nonrelativistic mass defined by 
Eq. f[3"3"j) parameterizes the momentum dependence of the single particle potential and is the 
result of a quadratic parameterization of the single particle spectrum. On the other hand, 
the relativistic Dirac mass is defined by Eq. (|30|) through the scalar part of the nucleon 
self-energy in the Dirac field equation which is absorbed into the effective mass. 

IV. SATURATION 

In this section the saturation mechanism and properties of nuclear matter in relativis- 
tic many-body approaches will be discussed. This saturation mechanism is quite different 
compared to nonrelativistic theories. In contrast to nonrelativistic theories, one has large 
attractive scalar fields and repulsive vector fields of several hundred MeV. Saturation occurs 
through the interplay between these large attractive and repulsive fields. The small nuclear 
binding energy arises from the cancellation between the large attractive scalar field and the 
large repulsive vector field. 

In RMF theory the large attractive scalar field Ssjs generated by the cx-meson, and the 



repulsive vector field E originates from the cu-meson |31|, |48j . In these relativistic approaches 
the vector field grows linearly with density whereas the scalar field saturates at large densities 
which leads finally to saturation. Correlation effects and Fock exchange terms are not 
present in these RMF models. However, this does not mean that the saturation mechanism 
is dominated by the mean-field or Hartree contribution and that exchange terms and higher 
correlations play only a minor role. In microscopic many-body approaches such as BHF and 
DBHF calculations based on realistic interactions, correlations and Fock exchange terms 
play an important role in the saturation mechanism. Therefore, RMF models which use the 
coupling strength values of the realistic interaction can qualitatively reproduce saturation, 
but completely fails in a quantitative description. Therefore, the parameters of RMF models 
have to be obtained from a fit to properties of finite nuclei and nuclear matter or are 
determined by a fit to the self-energies of a microscopic many-body calculation such as the 
DBHF approach. 

The binding energy results of some relativistic microscopic many-body approaches, a 
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DBHF calculation based on projections techniques 37J, |40|, |4l| and one based on the fit 



method 



391 ]. are depicted in Fig. [6) The saturation behavior of three different RMF theories 



are shown in Fig. EJ They are the NL3 [69] , DD [70] , and DDRMF [56] parameterizations. 
The parameters of the NL3 and the density dependent coupling functions of the DD model 
were obtained from a fit to properties of finite nuclei. The parameters of the DDRMF model 
are determined by a fit to the self-energies of a DBHF approach with a slight modification 
at saturation density to better reproduce properties of finite nuclei. In addition, nonrela- 
tivistic BHF calculations from Ref. [t^] based on the local Argonne Vis 78] are shown with 
and without three-body force. The included microscopic three-body force is based on a 
meson exchange model [^j]. The curves of all the microscopic and RMF methods are very 
similar below a density of approximately 0.2 fm -3 . This fact indicates that these are the 
density ranges which are at present reasonably well controlled by state-of-the-art many-body 
calculations. At higher densities noticeable differences occur. 

Another observation that becomes evident from fig. [6l by comparing microscopic calcu- 
lations to RMF calculations fitted to finite nuclei is that in nuclear matter the BHF and 
DBHF calculations lead to more binding than the RMF calculations. The prediction of a 
soft EoS is the general outcome of a microscopic many-body calculation. Recent Quantum 
Monte Carlo calculations for symmetric nuclear matter [79] show the same tendency. It 
should be noticed that this observation is supported by corresponding observables extracted 
from heavy ion reactions, where supranormal densities up to about three times saturation 
density are probed. Heavy ion data for tranverse flow [8(J or from kaon production 
support the picture of a soft EoS in symmetric nuclear matter. The compression moduli of 
the RMF models, which are presented in Tabled! are too high compared to the experimental 
value of K — 225 ± 15 MeV, except for the DD model where the nuclear incompressibility 
was fixed to K — 240 MeV. In general, RMF fits to finite nuclei req uire a relatively high 
compression modulus K of about 300 MeV at saturation density 481 ] . However, at higher 
densities the DDRMF model has the softest EoS, since this special procedure for the DD 
model can not avoid that at high density the EoS becomes very stiff. 

Furthermore, the models NL3 and DD yield rather low saturation densities around 0.15 
fm~ 3 whereas DDRMF saturates at a rather high density of p = 0.178 fm~ 3 . The saturation 
densities of both DBHF calculations are rather high, in particular for the calculation based 
on the fit method with a saturation density around 0.19 fm~ 3 , whereas the saturation binding 
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FIG. 6: Comparison of several EoSs from relativistic many-body approaches. The relativistic ab 
initio approaches are a DBHF using the fit method and one using projection techniques. The 
relativistic mean-field models presented are the nonlinear NL3 model, the DD model, and the 
DDRMF model. In addition, a nonrelativistic BHF is shown with three-body force and without 
three-body force. 

energies are in good agreement with empirical values. In addition, it becomes evident from 
Fig. [6] that, although nonrelativistic BHF calculations were able to describe the nuclear 
saturation mechanism qualitatively, they failed quantitatively compared to the relativistic 
DBHF calculations. Systematic studies for a large variety of NN interactions showed that 
saturation points were always allocated on a so-called Coester band in the E/A — p plane 
which does not meet the empirical region of saturation, which is demonstrated in Fig. [7J In 
particular for modern one boson exchange potentials, these nonrelativistic calculations lead 
to as well too much binding as too large saturation densities. The relativistic calculations 
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NL3 [69] 


DD [70] 


DDRMF [56J 


DBHF [39] DBHF [40, 41] 










(fit) 


fproi.) 


Po [f m 3 ] 


0.1482 


0.1487 


0.178 


« 0.19 


0.181 


E/A [MeV] 


-16.24 


-16.021 


-16.25 


-16.7 


-16.15 


K [MeV] 


271.5 


240 


337 


233 


230 


a s [MeV] 


37.4 


31.6 


32.1 


« 30 


34.36 



TABLE I: Saturation properties of nuclear matter are compared for the NL3 model, the DD 
model, the DDRMF model, a DBHF approach based on the fit method, and a DBHF approach 
based on projection techniques. The quantities listed include the saturation density po, the binding 
energy E/A at saturation density, the compressibility modulus K and the symmetry energy as at 
saturation density. 

do a much better job, which can be observed in Fig. [7] for DBHF results taken from the 
calculations of Brockman and Machleidt 



improved techniques from Ref. 



29| and from the more recent calculations based on 
37j using the Bonn potentials. The inclusion of three-body 
forces in the nonrelativistic BHF calculations makes the EoS stiffer, since the contributions 
from these forces are in total repulsive. Therefore, the inclusion of the three-body forces in 
these nonrelativistic BHF calculations lead to a shift of the saturation point away from the 
nonrelativistic Coester band towards its relativistic counterpart where the DBHF results are 
allocated. This shift of the saturation point is demonstrated for the BHF calculations [l?! 
denoted by triangles in Fig. where saturation density is shifted from a value of 0.265 fm -3 
to a value of 0.198 fm -3 . 

One may get the impression that the three-body forces, which are required in nonrelativis- 
tic calculations to achieve the empirical saturation point, corresponds to a parameterisation 
of the relativistic features, which are included in DBHF calculations. Indeed, if one ex- 
pands the propagator of a nucleon in the nuclear medium in terms of the propagators for 
the free nucleon, the lowest order term, which represents a change of the Dirac mass or 
the enhancement of the small component in the in-medium Dirac spinor, which is described 
by this in medium mass, is displayed in Fig. |HJ This so-called Z-graph involves two extra 
nucleon in addition to the propagating one and therefore should be represented in a nonrel- 
ativistic theory, which ignores this medium dependence of the Dirac spinor, by a three-body 
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FIG. 7: Saturation points of relativistic DBHF calculations (full symbols) are compared to several 
nonrelativistic BHF calculations without three-body forces (open symbols) and to the nonrela- 
tivistic BHF calculation including three-body forces (dotted triangular). The saturation points of 
the relativistic DBHF calculations are from Brockman and Machleidt (BM) and of the Tuebingen 
group. 
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The energy functional of nuclear matter can be expanded in terms of the asymmetry 
parameter (3 = (p n — p p ) / p (p n and p p are the neutron and proton densities, respectively) 
which leads to a parabolic dependence on 



E b {p, 13) = E{p) + E sym (p)f3 2 + 0(P 4 



(37) 



The first quantity, E(p), yields the binding in isospin symmetric nuclear matter, whereas 
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FIG. 8: Modification of the propagator representing the enhancement of the small component of 
the nucleon Dirac spinor in the medium due to the interaction with other nucleons. 

the quantity E sym (p) is called symmetry energy and characterizes the isospin dependence of 
the nuclear equation of state (EoS). Therefore, it is denned as 



E sym (p) = - 



d 2 E b (p,(3) 



(38) 

0=0 



d/3 2 

The density dependence of this symmetry energy for the models presented in Table [J is 
shown in Fig. [9j It can be observed that the symmetry energies in the various models are 
rather similar at a density near 0.1 fm -3 . In phenomenological models, such as RMF mod- 
els, the symmetry energy is constrained by the skin thickness of heavy nuclei which seems 
to fix the symmetry energy at an average density of about 0.1 fm -3 . The results of the 
microscopic approaches coincide at this density with RMF phenomenology demonstrating 
that the low density behavior of the microscopic calculation is in agreement with the con- 
straints from finite nuclei. However, at densities above the saturation density the predictions 
for the symmetry energy are quite different for the various many-body calculations. This 
high density behavior of the symmetry energy is essential for the description of the struc- 
ture and of the stability of neutron stars. Although microscopic DBHF calculations give 
a relatively soft EoS in isospin symmetric nuclear matter, these calculations yield a rather 
stiff isospin dependence, respectively symmetry energy, at high density. In general, they 
are significantly stiffer than in non-relativistic BHF approaches, such as the extended BHF 
(EBHF) approach Tlj] using the Argonne V14 jsj. On the phenomenological side, the NL3 
model has the stiff est EoS in Fig. |9] and the symmetry energy rises almost linearly with 
the density. In contrast, the DD model exhibits a considerable flattening. Therefore, the 
different phenomenological relativistic mean-field models can yield as well a soft as a stiff 
isopin dependence of the EoS. 
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FIG. 9: Comparison of the symmetry energy obtained from relativistic many-body calculations. 
The relativistic ab initio approaches are a DBHF using the fit method and one using projection 
techniques. The relativistic mean-field models presented are the nonlinear NL3 model, the DD 
model, and the DDRMF model. In addition, the result of a nonrelativistic extended BHF is 
shown. 



V. FINITE NUCLEI 



Many investigations have been devoted to a relativistic description of nuclei. Most of 
these studies, however, have been performed employing phenomenological models for the 
NN interaction. A very brief summary on such studies will be presented in Sec. IV Al Until 
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now no relatvistic calculation for finite nuclei has been presented, which is based on a realistic 
NN interactions and evaluates both, the correlation effects and the Dirac effects directly 
for the finite system under consideration (see e.g. the review by N. Van Giai et al. 84]). In 
Sec. IV Bl we are going to present the status of the approximation schemes for such DBHF 
calculations in finite nuclei. Special attention will be paid to the spin orbit splitting of the 
single-particle spectrum in Sec. IV C\ as the spin-orbit term can be considered as a kind of 
fingerprint for relativistic effects in nuclear systems. 



A. Phenomenological Methods 



Very popular phenomeno 



as the RMF theory 
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ogical models are relativistic density functional theories, such 
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53 



69 



85j. This RMF framework, or in other words, the 



relativistic Hartree approach with the no-sea approximation is already discussed in Sec. Ill Bl 
This approach, however, is missing the pion exchange processes due to the absence of Fock 
terms. This missing pion plays an essential role in determining the isospin dependence 
of the shell evolutions [86]. Therefore, considerable effort has been devoted to the rel- 
ativistic Hartree-Fock (RHF) theory to include the Fock terms 
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87H91|. However, it 



failed in a quantitative description of these nuclei for a long time, because of its numer- 



ical comp 
functions 



exity. After the introduction of microscopic based density dependent coupling 
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94j . phenomenological density dependent coupling functions were recently in- 
troduced in the relativistic Hartree (RH) theory [5^] and in the relativistic Hartree-Fock 



(DDRHF) theory [95|, |96j. It has lead to a good quantitative description of some nuclei, as 



can be observed in Table [TTJ The good agreement is not so surprising, since the parameter- 



ization PKAI 



961 ] was fitted to these nuclei. It is a general feature that the parameters in 



these purely phenomenological theories have been adjusted to describe the saturation prop- 
erties of nuclear matter or the properties of stable nuclei located in the valley of /3 stability. 
Therefore, the predictive power of these phenomenological interactions is rather limited, in 
particular for exotic nuclear systems such as the neutron star crust and nuclei far away from 
the line of /3 stability. 



30 







PKA1 [96] 




exp. 


[97-1001 






(a, us, it, p) 










E/A [MeV] 




r c [fm] 


E/A [MeV] 


r c [fm] 


WQ 


-7.94 




2.80 


-7.98 


2.74 


40 Ca 


-8.54 




3.53 


-8.55 


3.48 


48 Ca 


-8.67 




3.49 


-8.67 


3.47 


90 Zr 


-8.71 




4.28 


-8.71 


4.27 


208 pfo 


-7.87 




5.51 


-7.87 


5.50 



TABLE II: Results for closed shell nuclei applying the purely phenomenological DDRHF parame- 
terizations PKA1. Experimental values are also given. 

B. Microscopic Based Methods. 

Ab initio approaches, such as the DBHF approach, are based on high precision realis- 
tic free space nucleon-nucleon interactions and the nuclear many-body problem is solved 
without any adjustment of parameters. Therefore, if calculations of this kind are successful 
in describing properties of normal nuclei, we have got a tool which connects properties of 
baryons in the vacuum (NN scattering data) with nuclear systems at densities around the 
saturation density of nuclear matter. Such a tool should have a high predictive power, when 
it is used in extreme cases such as nuclear matter at densities beyond the saturation point 
or in a highly isospin asymmetric nuclear environment. Although DBHF calculations have 
been quite successful in describing nuclear matter, and therfore it fulfills this requirement, 
full Dirac Brueckner calculations are still too complex to allow an application to finite nu- 
clei at present. Therefore, different approximation schemes have been developed, which will 
treat either the correlation effects or the relativistic effects in an approximative way. These 
approximation schemes shall be described in this section. 

1 . Treatment of the Correlation Effects in a Local Density Approximation. 

In the first approximation scheme, the Dirac effects are treated directly for the finite 
nucleus, whereas the correlation effects are treated in a local density approximation. This 
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treatment of the correlation effects can be accomplished by defining effective meson exchange 
interactions based upon the nuclear matter T matrix. The coupling constants of such an 
effective interaction, which are adjusted so as to reproduce the in-medium T matrix, are 



density dependent, since the in- medium T matrix is density dependent |92|, |93|. In this 
way, a semi-phenomenological relativistic density functional can be constructed in contrast 
to the purely phenomenological relativistic mean-field model. These semi-phenomenological 
relativistic density functionals can be divided into semi-phenomenological DDRH theories 
and semi-phenomenological DDRHF theories. 

In order to properly parameterize the DBHF self-energy components in isospin asym- 
metric nuclear matter, the coupling functions in the semi-phenomenological DDRH theory 
should be based on four different channels: scalar isoscalar, vector isoscalar, scalar isovec- 
tor, and vector isovector channel. These channels correspond to phenomenological exchange 



bosons, i.e. the a-, u-, 5-, and p-mesons and these coupling functions are given by 



with p s 
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y/mf + k 2 



(43) 



and 



(27T) 



d 3 k 



k 3 
3tt 2 



(44) 



are, respectively, the scalar and vector density of the particle (i = n,p) 



In initial DDRH calculations, such as in Refs. [58|, |94|, the isovector scalar 5 coupling is 
missing. However, this coupling provides a mechanism to account for the differences in the 
scalar self-energies and in the corresponding effective Dirac masses for neutrons and protons 
in isopin asymmetric nuclear matter 3^|. Therefore, it has important consequences for the 



dynamics of neutron-rich nuclei 



36|]. 
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FIG. 10: The isoscalar effective coupling functions g a and are plotted as a function of the baryon 
density for different values of the asymmetry parameter f3. 



The coupling constants of Eqs. (|39l)-(l42l) depend on the density p and on the proton- 



neutron asymmetry density p 3 . In practice, the dependence on /? 3 o 
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the coupling functions is 
5511 . The reason is 



40 



41 



in general ignored to keep the functional as simple as possible 
the weakness of this p 3 dependence. In Fig. [10] this weak dependence on p 3 can be observed 
for the isoscalar coupling functions, which are obtained from the self-energy components 



of the DBHF approach in Ref. 



41] . Furthermore, these coupling 



unctions decrease with 



increasing density, which is attributed to the correlation effects 93[.The results for the 
isovector coupling functions are presented in Fig. [TT1 

However, one observes clear deviations between the original DBHF EoS and the DDRH 
EoS based on this DBHF EoS via the parameterization ( 139]) - ([42]) . These deviations are not 
due to the dependence on p 3 , as is shown in Fig. [12] Therefore, the density dependent 
coupling functions should be extracted more carefully as has been done in the "naive" 
definitio n (j39l) -(T42"l) by a kind of renormalization. This fact has already been pointed out in 



Ref. 



41 



55]. The reason of these deviations are that two essential differences between the 
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FIG. 11: The isovector effective coupling functions g$ and g p are plottted as a function of the 
baryon density for different values of the asymmetry parameter fi. 



DBHF approach and the RMF theory concerning the structure of the self-energy, respectively 
the mean- field, exist. 

The first difference is the appearance of a spatial contribution of the vector self-energy Ey 
in the DBHF theory, which is not present in the RMF model. This self-energy component 
originates from Fock exchange contributions which are not present in the RMF theory. 
Therefore, this spatial Sy component has to be included in a proper way for an accurate 
reproduction of the properties of the original DBHF EoS. The effects of the Ey component in 
the Dirac equation for homogeneous nuclear matter can be absorbed into a renormalization 



of the scalar and time-like vector component of the self-energy [4l|, |56|. Firstly, the Sy 
component can be absorbed into the effective mass according to Eq. ([TO]) and this reduced 
effective mass has to be identified with the RMF effective mass, i.e. 



M + S s , i{k Fi ) =M E DDRMF_ (45) 
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FIG. 12: Comparison between results of the DBHF approach (solid line), the DDRMF approach 
(dashed line), and the renormalized DDRMF approach (dashed-dotted line) for pure neutron (up) 
and symmetric (down) matter. The DDRMF approach and the renormalized DDRMF approach 
presented are based on the DBHF approach using a spline interpolation. 



This leads to the renormalized scalar self-energy component 



In a corresponding way, however, using the energy density instead of the effective mass, the 
following expression for the normalized vector self-energy component is obtained 



st^ddrmf v /. \ ^v,i{k Fi )\SE Fi pi + m*p Stii 

^0,i - ^oA^Fi) 7^ • 



These renormalized self-energy components are now inserted into Eqs. f l39|) - f|42|) to obtain 
the renormalized density dependent coupling functions. This renormalized DDRMF theory 
yields nuclear matter properties, which are closer to those of the original EoS 41] as can be 
observed in Fig. [121 Such a theory also provides good results for nuclei, as can be seen in 
Table HTT1 for the renormalized DDRMF theory in Ref. 56]. 
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DDRMF [56J 






E/A [MeV] 


r c [fm] 


Wq 


-8.35 


2.67 


40 C a 


-8.73 


3.35 


48 C a 


-8.73 


3.39 


90 Zr 


-8.74 


4.12 


208 Pb 


-7.87 


5.30 



TABLE III: Results for closed shell nuclei applying the DDRMF parameterization from the DBHF 
results. 

The second difference is that the DBHF self-energy terms explicitly depend on the mo- 
mentum of the particle, a feature which is absent in RMF. This momentum dependence 
reflects the non-locality of the DBHF self-energy terms, which originates not only from the 
Fock exchan ge t erms but also from non-localities in the underlying NN interaction. In the 
work of Ref. 55| . a method was introduced to compensate for this momentum dependence. 
Therefore, the full DBHF self-energies are expanded around the Fermi momentum, 



E(Jfc, hp) = £(JfeF, fcp) + (k 2 - k F f^ kp) 



dk< 



+ 0(^ 4 ). 



(48) 



k=k^ 



The first term in Eq. (148]) can be regarded as the Hartree self-energy, whereas the second 
term is a measure of the momentum dependence around the Fermi surface. In the next 
step, one should identify the energy density of DBHF EoS with that of the DDRH theory. 
Furthermore, the momentum will be integrated out by averaging over the Fermi sphere, 

4 



P^ DRH {k ¥ ) 



(27t)- J\ k \<L b 

pE G (/c, k F 



d 6 kJ: o (k : k F ) 
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E Q (k F ) 



Instead of extracting 



k F ) 



dk 2 



(49) 



(50) 



k=k-p 



directly from DBHF calculations, Xq(^f) is calculated numerically by adjusting the DDRH 
binding energy to that of the DBHF EoS in Ref. 55J. Therefore, this correction of the self- 
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energy can be translated into a modification of the density dependent coupling functions, 



~g\k ¥ ) = g 2 (k F ) 



(51) 



3 ~ r S G (fc F ) 

Making the assumption that the ratio only depends weakly on the Fermi momentum, 

one can introduce the following momentum corrected coupling functions, 



9<x(fa) = 9a{k F )\h-- Cakl, 



(52) 



where the constants ( a are determined by fitting to the DBHF EoS. This momentum de- 
pendence correction improves the calculated results as can be observed from Table IIVI for 
some closed shell nuclei. 





DDRH [55J 




cor. DDRH [55] 






E/A [MeV] 


r c [fm] 


E/A [MeV] 


r c [fm] 


m Q 


-8.58 


2.75 


-8.30 


2.79 


40 C7a 


-9.02 


3.46 


-8.69 


3.50 


48 C7a 


-8.96 


3.49 


-8.63 


3.53 


90 Zr 


-8.98 


4.26 


-8.63 


4.31 


208 pfo 


-8.17 


5.53 


-7.82 


5.60 



TABLE IV: Results for closed shell nuclei applying the DDRH parameterizations without and 
with a correction for the momentum dependence. These parameterizations are based on DBHF 
calculations using Bonn A. 



These essential differences between the DBHF approach and the RMF theory concerning 
the structure of the self-energy mainly originate from the omission of the Fock exchange 
terms. An improvement, therefore, was the density dependent relativistic Hartree-Fock 
(DDRHF) theory with its inclusion of the Fock terms. Since Fock terms are present, one 
obtains a spatial contribution of the vector self-energy in the DDRHF theory and mo- 
mentum dependent self-energy components as in the DBHF approach. Up to now the rela- 
tivistic DBHF in-medium T matrix is presented by an effective interaction Lagrangian in the 
DDRHF theory, which only contains the a-,u-,7T-, and the p-meson {q2|. The corresponding 
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self-energy components in this DDRHF theory are 

1 rk F j t3 m * 2 

j <J>J 

m *2 _|_ m *2 \-2k*^Q* 
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q z (m* + m*) 2 

with the meson propagator, 

Afefo g) = -71 T77 x 9 • (56) 

v ' -(k*» - q*i*)(k* - q*) + rn 2 k 

A typical feature of this theory compared to Hartree versions is the presence of a pion 
exchange potential. This pion exchange potential does not appear in Hartree theories, since 
this potential only contributes to the Fock terms. Nevertheless, these DDRHF models based 
on microscopic approaches are not so successful yet as DDRH models, as can be concluded 
from Table [V] A possible reason is that only the isoscalar coupling functions are density 
dependent, whereas the coupling functions of t he n- meson and the p- meson remain density 



independent in these DDRHF theories 45l . Il0l| . Another observation is that the 5-meson is 



not present in these DDRHF theories. Therefore, these theories can be further improved by 
including the 5-meson and extent the density dependence to all coupling functions, which 
has not been done until now for DDRHF theories based on microscopic calculations such as 
the DBHF approach. 



38 



2. Treatment of the Dirac Effects in a Kind of Local Density Approximation. 



In this scheme 4J, Il02[ Il03j , the information about the structure of the Dirac spinors is 
taken from the investigations of nuclear matter. This means that the Dirac effects are taken 
into account via a kind of local density approximation (LDA). Therefore, this approximation 
scheme will be called Dirac LDA. On the other hand, the correlations effects are taken 
into account by solving the Bethe-Goldstone equation directly for the finite nucleus under 
consideration. This means that the self-consistency requirements of a conventional BHF 
calculation for finite nuclei are fulfilled, whereas relativistic effects are taken into account 
by evaluating the potential matrix elements in terms of in-medium Dirac spinors. 



In this Dirac LDA scheme 



44[ , the calculated binding energy and the radius of the charge 



distribution is increased compared to the results of nonrelativistic BHF calculations [44| . as 
can be observed in Fig. [13] in case of the isospin symmetric 16 O nucleus. Therefore, the 
results of these Dirac LDA calculations are generally in better agreement with experiment 
than the ones of the nonrelativistic BHF calculations, in which the medium dependence 
of the Dirac spinors is ignored. In addition one can see from Fig. [HI] that the results of 
this Dirac LDA approach are very close to those obtained within the DDRHF calculation. 
The same feature has also been observed for other isospin symmetric nuclei. This can be 
considered as a mutual support of these two approximation schemes indicating that LDA 
is a good approximation for treating correlation effects as well as Dirac effects. The result 





E/A [MeV] 


HF2 [45J 

(a, u), 7r) 


r c [fm] 


E/A [MeV] 


RDHF3B [1011 


r c [fm] 


16 


-7.73 




2.48 


-7.57 




2.59 


40 Ca 


-8.09 




3.14 


-7.94 




3.26 


48 Ca 


-7.90 




3.16 


-7.70 




3.27 


90 Zr 








-7.78 




3.99 


208p b 








-6.61 




5.17 



TABLE V: Results for closed shell nuclei applying the DDRHF parameterizations from DBHF 
results. 
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FIG. 13: The binding energy per nucleon and radius of the charge distribution obtained from a 
nonrelativistic BHF approach and some approaches based on the DBHF results are compared to 
the experimental value. These approaches based on DBHF results are the Dirac LDA, the density 
dependent relativistic Hartree (DDRH) theory, and the density dependent relativistic Hartree-Fock 
(DDRHF) theory. 



of a complete DBHF calculation should be in the vicinity of the Dirac LDA and DDRHF 
results. 

On the other hand it can be observed in Fig. [13] that the result of the semi- 
phenomenological DDRH calculation of Ref. 44] is in much better agreement with the 
experimental value than that of the Dirac LDA calculation. However, this good agreement 
with the experimental value seems just fortuitous, since it is lost when Fock terms are 
included. 

It should be recalled that all these calculations discussed in this subsection ignore rear- 
rangement terms, although e.g. the density dependence of the effective coupling constants 



40 



in the DDRHF approximation scheme would call for such rearrangement terms. They have 
been dropped because the aim of this study has been to develop a reliable approximation 
scheme for DBHF in finite nuclei. The inclusion of rearrangement terms or the correspond- 
ing extended microscopic definition of the nucleon self-energy would be one possible step 
for future investigations to improve the results of relativistic calculations beyond the DBHF 
approach. 



C. Spin Orbit Splitting 



In atomic physics, the spin orbit splitting of the single-particle spectrum is a genuine 
relativistic effect. Therefore it is quite natural to consider the spin orbit splitting in nuclear 
system as a possiblle candidate to probe the relevance of relativistic features for the nuclear 
many-body system. Therefore this section is devoted to the discussion of the spin orbit 
term. Reducing the Dirac equation for a nucleon in a RMF theory to a non-relativistic 
Schrodinger equation one obtains a spin orbit term of the form 

1 D'{r) s • 1 



U ls (r) 



2m D{r) r 



(57) 



with 



D(r) =m + S s (r) + £ + £ (r). (58) 

From the term s • 1 in Eq. (!57|) . it is expected that the spin orbit splitting tends to become 
larger with increasing orbital angular momentum. However, inspecting the experimental 
data for the example of the nucleus 16 O, one finds that the splitting of the p orbits is larger 
than the splitting of the d orbits as can be observed from Table IVIl 





BHF 


BHF + 2hlp + lh2p 


BHF + 2hlp + lh2p 


experiment 








+ Dirac effects 




Ae p 


3.95 


4.11 


6.19 


6.18 




5.80 


5.50 


5.50 


5.09 



TABLE VI: Spin orbit splitting in 16 O. The spin orbit splitting for the d shell and the p shell in 
the various approximation schemes are given in MeV. 
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Nonrelativistic BHF calculations can not reproduce these experimental results for the 
spin or bit sp litting, i.e. it predicts a larger spin orbit splitting for the d shell than for the 
p shell 104l | . It is well known that contributions beyond the BHF definition of the single- 
particle potential (rearrangement terms) can modify the spin orbit splitting to quite some 
extent. Therefore the BHF contribution has been supplemented with the 2-hole 1-particle 
(2hlp) term and with the 2-particle 1-hole (2plh) term. The effect of the 2hlp contribution is 
repulsive and essentially cancels the effects of the 2plh term, which is attractive. Ther efore , 
the inclusion of these terms do not change the spin orbit picture in a significant way 104 1. 

The Dirac LDA method is used to investigate the influence of the relativistic effects on the 
spin orbit splitting. The inclusion of these relativistic effects lead to a drastic enhancement 
of the spin orbit splitting for the p she ll, whereas the Dirac effects do not really enhance the 
spin-orbit splitting of the d shell 104J, as can be deduced from Table IVT1 



A simple explanation for these findings c an be pro vided: The spin orbit splitting is to the 
first order proportional to the inverse of m* 105l . Il06| . since the value of D(r) in Eq. ([58]) can 
be approximated by 2m* (r). In Fig. [3] we have displayed the density dependence of the Dirac 
mass. So its evident that this Dirac mass, m*, is much smaller in the dense interior of the 
nucleus than at the lower densities at the surface. Therefore, the spin orbit splitting of the 
p shell is much stronger enhanced than the one of the d shell, the latter having the largest 
amplitudes at the low density region near the surface. Thus, this calculation including the 
relativistic effects predicts that the spin orbit splitting in 16 for the p shell is larger than 
the one for the d shell, which is in agreement with experimental measurements. 

Of course there also exist other approaches which reproduce the spin-orbit effects. As 
an example we men tion the non-relativistic many-body calculation on 16 by Pieper and 



Pandharipande 



1071 ] . Also they can reproduce the spin-orbit splitting of the p-shell, if the 



effects of a 3-body force are taken into account. Again we have this feature that effects 
obtained within the Dirac phenomenology are described within a non-relativistic approach 
using a 3-nucleon force. It is not clear, however, whether this 3-body force will also lead to 
a proper result for the splitting in the <i-shell. 
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VI. SATURATION WITH LOW-MOMENTUM INTERACTIONS AND RELA- 
TIVISTIC EFFECTS. 



As we have already discussed in the preceeding parts it is still a challenge for theoretical 
nuclear physics to predict the properties of finite nuclei from the knowledge of a bare NN 
interaction, which is fitted to describe the two-nucleon data. The DBHF approach, which 
yields a rather good description of the saturation properties of symmetric nuclear matter 
without the necessity to introduce a three-nucleon force, has not yet been extended to be 
applicable also to finite nuclei. The reason for this missing extension is the fact that it is 
not easy to evaluate the correlation effects as well as the relativistic effects consistently for 
finite nuclei (see Sec. IV Bl and Ref. |84j). 

A possible way out of this problem is the use of realistic interactions, which are renor- 
malised to be used in a Hilbert space with low relative momenta between the interacting 
nucleons. In this way all short-range components of the NN intteraction are integrated out 
and there is no need for a non-perturbative treatment of NN correlations. Such r estrictions 



of the nuclear structure calculation to the low- momentum components 22J, 123, Il08l . Il09 | 
have become very popular during the last year within the nonrelativistic framework. In 
this section we would like to present first studies, which employ these technique also for the 
relativistic approach. 

A. The Low-Momentum Interaction. 

The idea of restricting the nuclear structure calculation to the low-momentum compo- 
nents is very reasonable. The long-range or low-momentum part of the NN interaction 
is fairly well described in terms of meson-exchange, whereas the quark degrees of freedom 
are more important in describing the short-rang or high-momentum components of the NN 
interaction. Therefore, it is quite attractive to disentangle these low-momentum and high- 
momentum components from each other to construct a low-momentum potential. This 
means that a model space should be defined, which accounts for the low-momentum de- 
grees of freedom and renormalizes the effective Hamiltonian for this low-momentum regime 
to account for the effects of the high-momentum parts. One of the most used methods to 



disentangle these parts is the unitary-model-operator approach (UMOA) [26]. 
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Therefore, one has to define two projection operators. As in all well-known model space 
techniques, the operator P projects onto the model space of two-nucleon wave functions 
with relative momenta k smaller than a chosen cut-off A, whereas the operator Q projects 
onto the complement of this subspace, the high-momentum subspace. In addition, these 
operators P and Q satisfy the usual relations like P + Q = 1, P 2 = P, Q 2 = Q, and 
PQ = = QP. In the UMOA approach, the aim is to define a unitary transformation U in 
such a way that the transformed Hamiltonian does not couple P and Q, i.e., 

QU^HUP = 0. (59) 



The technique to determine this unitary transformation has nicely been described in Ref 27 



110 1 . With the help of this unitary transformation an effective Hamiltonian H e ff = ho + V e ff 
can be obtained. It contains the kinetic energy /z and an effective interaction V e ff given by 

V eff = V lowk = U-\ho + V)U - ho. (60) 

The eigenvalues, which are obtained in the diagonalization of this effective Hamiltonian 
in the P-space, are identical to those, which are obtained by diagonalizing the original 
Hamiltonian H = ho + V in the complete space. This means that solving the Lipmann- 
Schwinger equation for NN scattering using this V e ff with a cut-off A yields the same phase 
shifts as obtained for the original realistic interaction V without a cut-off. Therefore, this 
V e ff, which from now on will be called Vi ow k, can also be regarded as a realistic interaction 
like e.g. the CD Bonn or Argonne V18 interactions, since it reproduces the NN scattering 
phase shifts. 

An interesting feature is that the resulting V[ ow k is found to be essentially independent of 
the underlying realistic interaction V, when V is fitted to the experimental phase shifts and 
the cut-off A is chosen around A = 2 fm _1 . Hence, one is able to obtain a low- momentum 
potential V[ ow k, which in a model independent manner describes the low- momentum com- 
ponent of realistic NN interactions in a more or less unique way. 

Furthermore, the Vi ow k does not induce any short-range correlations into the nuclear 
wave function, because the high-momentum components have been integrated out by means 
of the unitary transformation of Eq. (1591) . Therefore, mean-field calculations using Vi ow k 
already lead to reasonable results and corrections of many-body theories beyond mean-field 



are weak 



27] 
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FIG. 14: Binding energy per nucleon of symmetric nuclear matter of a HF calculation (solid 
line) and of a BHF calculation (dashed dotted line) employing a V[ ow k interaction are plotted. In 
addition, the binding energy per nucleon of a nonrelativistic BHF with three-body force calculations 
(dashed double dotted line) and without three-body force (dashed line) is presented using the 
Argonne Vis potential. . 



B. Saturation Behavior of Vi ow k 

An important feature of an EoS is the saturation behavior. However, the saturation 
behavior of Vi ow k is problematic, since in HF calculations of nuclear matte r on e obtains a 



binding energy per nucleon increasing with density in a monotonic way |42l . 11 11] . as can be 
seen in Fig. [TH Hence, HF calculations employing Vi ow k do not exhibit a minimum in the 
energy as a function of the density and the emergence of a saturation point is prevented in 
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isospin symmetric nuclear matter 



27 



28 



This absence of saturation also leads 
to problems in the description of finite nuclei. In fact, Hartree-Fock calculations of finit e 
nuclei using Vi ow k yield higly compressed nuclei, if the variational space is not restricted [ill ]. 

This absence of saturation is one of the major problems of calculations employing Vi ow k. 
Even the inclusion of correlations beyond the HF approximation, e.g. by means of the BHF 
approximation 42] , can not solve this problem as can be observed from Fig. [TH Employing a 



conventional model for a realistic interaction, like e.g. the Argonne V18, which is not reduced 
to its low-momentum components, one obtains a saturation point in calculations which 



account for correlations beyond the mean- field approximation [77J (see Fig. [T4"j) . Therefore, 
one may argue that the lack of short-range correlation effects, which have been integrated out 
by means of the unitary transformation of Eq. floTJl) . prev ents the emergence of a saturation 



point in calculations of symmetric nuclear matter 28 



108j . This means that the Vi 



owk 



approach can not reproduce the saturation of nuclear matter as it misses the quenching of 
short-range correlations in the nuclear medium. This quenching of short-range correlations 
is included in sophisticated calculations employing one of the conventional models for a 
realistic NN interaction. 

Although calculations, like Brueckner-Hartree-Fock (BHF), for these conventional realis- 
tic interactions lead to a saturation point, they are not able to reproduce the experimental 
data as already mentioned in Sec. IIVI In order to o btai n a saturation point closer to the 
empirical data one has to include a three-body force jl!2 | or inclu de relatiy istic effects, e.g. 



within a Dirac-Brueckner-Hartree-Fock (DBHF) approach 40 



41 



103 



113j . Since it can be 



expected from these findings that these measures will also improve the satur ation behavior of 
calculations employing V[ ow k, this issue has been investigated. In Ref. [1081 ]. it is shown that 
the addition of three-body forces can lead to saturation in nuclear matter. The inclusion of 
the relativistic effects will now be elaborately discussed in Sec. IVI CI 



C. Relativistic Effects 



Relativistic effects lead to a successful microscopic description of the saturation properties 
of nuclear matter in the DBHF approach. The main reason for this success can be attributed 
to the fact that the matrix elements of the bare nucleon-nucleon interaction become density 
dependent, since the spinor contains the reduced effective mass [3^]. This density dependence 
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FIG. 15: Binding energy per nucleon of isospin symmetric nuclear matter of a HF calculation 
(dashed line) and of a BHF calculation (dotted line) employing a density dependent Vi ow k inter- 
action are plotted. In addition, a HF calculation using a standard Vi ow k interaction (solid line) is 
added to the figure. 



is absent in nonrelativistic Hartree-Fock or Brueckner calculations. 

These relativistic effects are introduced in Vi ow k leading to a density dependent 
Viowk{p) |42|. This is achieved by calculating for each density the underlying realistic in- 
teraction in terms of these dressed Dirac spinors. The density dependent Vi ow k{p) can then 
be obtained using the standard techniques given in Sec. lyTAl Results for the energy per 
nucleon of isospin symmetric nuclear matter as a function of the density p obtained from HF 
and BHF calculations employing density dependent V/ ow fc(p) are displayed in Fig. [T5j 
The density dependent Vi ow k{p) presented in this figure is based on the matrix elements of 
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the Bonn A potential using the Dirac spinors, which are appropriate for this density. There- 
fore, the medium properties of the nucleons, which are used to dress these Dirac spinors, are 



taken from the EoS presented in Refs. [40. l4l| . Employing this density dependent V[ ow k{p) 



interaction with a cut-off A of 2 fm" 1 , one already obtains a saturation point in the HF 
calculation 42^ . as can be observed from Fig. [T5J The reason is that the attractive contri- 
butions of the scalar meson cr-meson are reduced, because the small component in the Dirac 
spinor is enlarged due to the fact that the effective Dirac mass is smaller than the bare mass. 
Because this effect grows with increasing density, the HF calculation employing the density 
dependent Vi ow k{p) interaction already leads to a minimum in the energy as a function of 
the density. Although this HF calculation yields a saturation density in the neighbourhood 
of the empirical region, it still provides too little binding, i.e. about -12 MeV per nucleon. 

This HF result can be improved by taking into account the effects of NN correlations 
within the BHF approximation, which yields about 2 MeV more binding in the region of 
saturation due to the additional correlations from NN states. As in the case of the standard 
Vi ow k, these NN correlations beyond mean-field also are raher weak for the density depen- 
dent Vi ow k{p) compared to the conventional realistic interactions and they even completely 
vanish at high densities. It can be explained by the fact that only intermediate two-particle 
states with momenta below the cut-off A can be taken into account in the Bethe-Goldstone 
equation, which in the two-particle center of mass frame takes for an effective interaction 
with cut-off A the form 

A 

2. 



G(k', k, e k ) eff = V eff (k', k) + / q 2 dqV eff (k', q) 

Qp 







-G(q,k,e k ), (61) 

2e k - 2e q + vq 

where e« with i = k,q are single-particle energies and Qp is the Pauli operator. Since one has 
to take into account that V e ff is designed for a model space with relative momenta smaller 
than A, the integral in Eq. (l6Tj) is restricted to momenta q below the cut-off A. Therefore, 
the NN correlations become weaker at higher densities because of the lack of phase space 
for these correlations in contrast to the case of conventional realistic NN interactions. In 
the latter case, the upper integration limit is absent in Eq. f lBTT) . 

Another intriguing feature of the Vi ow k approach has been the finding that the resulting 
low-momentum interaction is rather insensitive to the NN interaction on which it is based. 
One obtains essentially a unique interaction model for the low-momentum regime. The same 

48 



30 



20- 



> 



10- 



< 

S3 Oh 



-10 r. 



-20 



"i 1 r 



"i 1 r 



V lowk (p) in BHF approximation 



Bonn A 
BonnB 



/ 

* 

* 

* 



♦ * 
♦ ♦ 



j i_ 



J L 



0.1 



0.2 



0.3 



-3-, 



Density [fm ] 



FIG. 16: Binding energy of isospin symmetric nuclear matter as a function of density from BHF 
calculations employing the Bonn A and Bonn B potentials for the density dependent Vi ow k(p) 
interaction. 



is also true using the Vi ow k approach within the relativistic framework as can be seen from 
Fig. [TH The energy versus density curve which are obtained in BHF calculations using 
density-dependent Vi ow k(p) interactions based on different meson exchange potentials are 
rather close to each other. The differences in the calculated energies are significantly smaller 
as those obtained when the underlying Bonn A and Bonn B potentials are used directly 
without the renormalization to Vi ow k- The differences seem to be slightly stronger than 
the corresponding model dependence in the standard Vi ow k interactions using the modern 
potentials, which fit the A^ A" data with high precision. However, one must keep in mind, 
that the potentials Bonn A and B did not fit the phase-shifts that accurately and indeed 
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a similar model dependence is also obtained applying the renormalization scheme to these 
potentials in the conventional way ignoring the relativistic effects. 

In short, HF calculations using a density dependent Vi ow k(p) can already lead to reason- 
able results and the NN correlations beyond mean-field are rather weak. These properties 
of density dependent Vi ow k(p) give a promising prospect for studying finite nuclei in cal- 
culations which are based on a realistic iViV interaction, treating correlation effects in a 
perturbative manner. This should be feasible also for the description of finite nuclei. 

VII. SUMMARY AND OUTLOOK 

We have discussed the status of relativistic nuclear many-body calculations in nuclear 
matter and nuclei. These kind of calculations are based on DBHF approaches or relativistic 
density functional theories such as the RMF theory. The relativistic density functional theory 
is a phenomenological approach. The parameters of these phenomenological models have 
been adjusted to describe properties of isospin symmetric nuclear matter and of nuclei in the 
valley of j3 stability. Therefore, a problem may exist concerning the reliability for systems 
at large densities and proton-neutron asymmetries. In contrast, the DBHF approach is an 
ab initio approach and the nuclear many-body problem is treated microscopically. Hence, 
these DBHF approaches have a high predictive power, since the predictions for the nuclear 
EoS are essentially parameter free. A typical characteristic of both these relativistic nuclear 
many-body approaches are the large scalar and vector self-energy components with a size 
of several hundred MeV in contrast to nonrelativistic approaches. These large components 
cancel each other to a large extent in determining the energy of the nucleons. The large 
attractive scalar component of the self-energy, however, yields a significant enhancement of 
the small Dirac components in the nuclear medium. 

This genuine relativistic effect is described in terms of the Dirac mass, which can only be 
determined from relativistic approaches. It is defined through the scalar part of the nucleon 
self-energy. This Dirac mass decreases continuously with increasing density. The DBHF 
calculations based on projection techniques predict a Dirac mass splitting of m* Dn < m* D in 
neutron-rich matter in contrast to DBHF calculations based on the fit method. However,the 
explicit momentum dependence has already been averaged out in the fit method, which 
makes this procedure ambiguous in isospin asymmetric nuclear matter. Furthermore, RMF 
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theories with the isovector 5-meson included predict a Dirac mass splitting of m* D n < m* D p . 

However, the expression of an effective nucleon mass has also been used within the context 
of nonrelativistic calculations. This nonrelativistic mass is a measure of the nonlocality of 
the single particle potential. These nonlocalities of the single particle potential can be due 
to nonlocalities in space or in time. The nonrelativistic mass shows a typical peak structure 
around k-p, which is due the energy dependence of the nucleon self-energy or its nonlocality 
in time. This peak structure is observed in relativistic as well as nonrelativistic calculations. 
As a consequence of this energy dependence the nonrelativistic mass is not a smooth function 
of the momentum in contrast to the Dirac mass. 

The Landau mass, i.e. the nonrelativistic mass evaluated at k — kp, initially decreases 
with increasing density like the Dirac mass. However, it starts to rise again at high values of 
the density. The nonrelativistic mass derived from the DBHF approach based on projection 
techniques in general shows a nonrelativistic mass splitting of m* NRn > m* NRp , which is 
opposite to its Dirac mass splitting of m* Dn < m* Dp . However, it is in agreement with 
results from nonrelativistic BHF calculations. The opposite behavior between the Dirac 
mass splitting and the nonrelativistic mass splitting is not surprising, since these masses are 
based on completely different physical concepts. 

An important bench mark for nuclear structure calculations is its success in describing the 
saturation in isospin symmetric nuclear matter. Although nonrelativistic BHF calculations 
were able to describe the nuclear saturation mechanism qualitatively, they failed quantita- 
tively. Systematic studies for a large variety of NN interactions showed that saturation 
points were always allocated on a so-called Coester band in the E/A — p plane which does 
not meet the empirical region of saturation. 

The relativistic DBHF calculations do a much better job. The saturation mechanism 
of these relativistic calculations is quite different compared to nonrelativistic theories. The 
density dependence of Dirac spinors described in terms of an effective Dirac mass decreasing 
with density leads to a density dependence of the meson exchange interaction terms, which 
determines the minimum in the energy versus density curve for nuclear matter. 

In nonrelativistic calculations the saturation point is shifted away from the Coester band 
towards the values of its relativistic counterpart only after the inclusion of three-body forces. 
Therefore, it is often argued that three-body forces in nonrelativistic calculations are required 
to simulate the dressing of the two-body interaction by in-medium spinors in the Dirac 
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phenomenology of the relativistic approaches. It should be noted, however, that there are 
other sources, beside these Dirac effects, which lead to three-nucleon forces. Examples are 
sub-nucleonic degrees of freedom like the medium modification of the meson propagators or 
the inclusion of excitation modes for the interacting nucleons. 

Another issue is the symmetry energy in relativistic microscopic nuclear many-body ap- 
proaches. Although the symmetry energy in nonrelativistic and relativistic models are rather 
similar at a density near 0.1 fm~ 3 , the predictions are quite different at higher densities. 
However, this high density behavior of the symmetry energy is essential for the description 
of the structure and of the stability of neutron stars. Microscopic DBHF calculations yield 
a rather stiff isospin dependence, respectively symmetry energy, at high density, whereas 
non-relativistic BHF approaches give a rather soft isospin dependence. 

Many investigations have been devoted to a relativistic description of finite nuclei. Very 
popular phenomenological models are relativistic density functional theories. It is a general 
feature that the parameters in these purely phenomenological theories have been adjusted 
to describe the saturation properties of nuclear matter or the properties of stable nuclei 
located in the valley of ft stability. Modern relativistic density functionals are able to yield 
a good quantitative description of these nuclear systems. However, the predictive power of 
these phenomenological interactions may be rather limited for exotic nuclear systems such 
as the neutron star crust and nuclei far away from the line of ft stability. 

Ab initio approaches such as the DBHF approach, which are based on high precision re- 
alistic NN interactions, are more ambitious. However, full Dirac Brueckner calculations are 
still too complex to allow an application to finite nuclei at present. Therefore, two different 
approximation schemes have been developed: Dirac LDA and the semi-phenomenological 
relativistic density functional theory. 

In the Dirac LDA approach, the Dirac effects are taken into account via a kind of local 
density approximation, whereas the correlations effects are taken into account by solving 
the Bethe-Goldstone equation directly in the finite nucleus under consideration. The results 
of these LDA Dirac calculations are generally in better agreement with experiment than 
the ones of the nonrelativistic BHF calculations, in which the medium dependence of the 
Dirac spinors is ignored. The relativistic effects also improve the microscopic understanding 
of details in nuclear structure like the spin orbit term. Further improvements may be 
achieved by including terms in the expansion of the nucleon self-energy, which go beyond 
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the BHF approximation. Such rearrangement terms are also required to lead to a symmetry- 
conserving approach. 

In the semi-phenomenological relativistic density functional theory, the Dirac effects are 
treated directly for the finite nucleus, whereas the correlation effects are treated in a local 
density approximation. The results are in better agreement with the experimental value 
than those from nonrelativistic BHF approaches and are comparable or even better than 
those from Dirac LDA. 

An alternative to these methods is to restrict the nuclear structure calculation to the 
low-momentum components of realistic AW interactions. An interesting feature of such a 
low momentum potential Vi ow k is that it is essentially independent of the underlying realistic 
interaction V, when V is fitted to the experimental phase shifts and the cut-off A is chosen 
around A = 2 fm _1 . However, within the nonrelativistic framework the the Vi ow k approach 
fails the test of predicting nuclear matter saturation. HF calculations employing Vi ow k do not 
lead to a saturation point in symmetric nuclear matter. This absence of saturation is one of 
the major problems of calculations employing Vi ow k, since even the inclusion of correlations 
beyond the HF approximation, e.g. by means of the BHF approximation, can not cure this 
problem. 

The inclusion of relativistic effects to obtain a density dependent Vi ow k(p) seems to be 
a promising prospect to solve this saturation problem, since HF calculations employing a 
density dependent Vi ow k(p) can already reproduce the saturation point in nuclear matter. 
Furthermore, this Vi ow k(p) does not induce any short-range correlations into the nuclear wave 
function, because the high- momentum components have been integrated out. Mean- field 
calculations using Vi ow k{p) already lead to reasonable results and corrections of many-body 
theories beyond mean-field like in the BHF approach are weak and can be treated in a 
perturbative way. Therefore, self-consistent relativistic HF calculations applying a density 
dependent Vi ow k(p) seem to be the future for studying finite nuclei in calculations which are 
based on a realistic NN interaction. 
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